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We define thermodynamic configurations and identify two primitives of discrete quantum processes 
between configurations for which heat and work can be defined in a natural way. This allows us 
to uncover a general second law for any discrete trajectory that consists of a sequence of these 
primitives, linking both equilibrium and non-equilibrium configurations. Moreover, in the limit of 
a discrete trajectory that passes through an infinite number of configurations, i.e. in the reversible 
limit, we recover the saturation of the second law. Finally, we show that for a discrete Carnot cycle 
operating between four configurations one recovers Carnot 's thermal efficiency. 



I. INTRODUCTION 

The intuitive meaning of heat and work in thermo- 
dynamics is that of two types of energetic resources, 
one fully controllable and useful, the other uncontrolled 
and wasteful. An impressive effort has been devoted 
to provide a consistent mathematical characterisation of 
these notions within a quantum mechanical description 
of physics PHI]- This is a challenge since in contrast to 
other thermodynamic quantities, such as internal energy 
and entropy, heat and work are not properties of indi- 
vidual states of a system. They are defined for continu- 
ous processes connecting different states [2, 8, 9 , imply- 
ing that their statistical fluctuations cannot be described 
in terms of a single system observable. Two-point cor- 
relation functions characterising the correlations along 
process paths are required - a problematic territory for 
quantum mechanics where definite trajectories cannot be 
fixed unless the system is continuously measured. Re- 
solving these issues has been the topic of a number of 
publications that have formulated quantum trajectory 
approaches [1 H El [SHIT] . 

In contrast here we focus on the mean values of heat 
and work where the analysis simplifies but still requires 
careful thought. We will adopt the identification of the 
system's internal energy with U(p) = tr[pH] where p 
is the density matrix describing the state of the sys- 
tem at given time, and H is its instantaneous Hamil- 
tonian. Clearly, a proper definition of this Hamiltonian 
is in general problematic! If the system is coupled to an 
environment the non-equilibrium behaviour of a general 
open system makes the definition of the system's Hamil- 
tonian ambiguous [TOl IT2TfT5] both, mathematically and 
experimentally. Ultimately the choice of the Hamiltonian 
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one assigns to the system must rely on the set of opera- 
tions and observables one can access experimentally. In 
many situations of physical interest H can be identified 
with the bare system Hamiltonian or an effective system 
Hamiltonian that incorporates the effect of the environ- 
ment. 

While the environment degrees arc in principle uncon- 
trolled, full control can be exerted over the temporal 
"variation" of the system Hamiltonian. For instance, the 
size of a container in which steam is pumped can be freely 
chosen and a piston can be attached to the container that 
can push the wheels of a train. A formal definition of 
mean heat and work is then obtained by considering an 
infinitesimal change of the internal energy 

dU = dtr[pH] = tr[pdH] + tr[dpH] , (1) 

associated with the time evolution of the system which 
brings its density matrix from p to p+dp while the Hamil- 
tonian varies from H to H + dH. The origin for the 
change of p may here be due to both, the variation of H 
induced by the experimenter and by the dynamics due 
to the coupling with the environment. The possibility 
of externally controlling H suggests to identify the first 
term on the r.h.s. of Eq. ([I]) with the average work done 
by the experimenter during the evolution. The second 
term describes the internal energy change due to a re- 
configuration of the system, i.e. a variation of the sys- 
tem's density matrix. This is an energy contribution over 
which the experimenter has no direct control and this is 
why it is associated with heat. The infinitesimal average 
heat absorbed by the system and the infinitesimal average 
work done on the system Ht)H2"5] are therefore defined 
as 

SQ := tv[dpH] and SW := tv[pdH] , (2) 

with the symbol 5 indicating that heat and work are in 
general no full differentials, i.e. they do not correspond 
to observables. 
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While the first law of thermodynamics states that the 
sum of the two average energy types is the average inter- 
nal energy, 

dU = 5Q + SW , (3) 

the split into these two types of energies is crucial as it 
allows the formulation of the second law of thermody- 
namics. A fundamental law of physics, it sets limits on 
the work extraction of heat engines and establishes the 
notion of irreversibility in physics. The second law can 
be phrased in form of Clausius' inequality: 

TdS > 5Q, (4) 

stating that the change in a system's entropy must be 
larger than the average heat absorbed by the system dur- 
ing a process. While the first law of thermodynamics is 
fundamental for any process, the second law was origi- 
nally stated for processes that start and end in equilib- 
rium. Recently, the non-equilibrium work relation due 
to Jarzynski has been used to argue that the second law 
should also hold for any process starting from equilib- 
rium, at temperature T, but ending in an arbitrary non- 
equilibrium state [21] • However, no conclusive argument 
has yet established the most general set of dynamical 
processes that obey the Clausius inequality [2"5] . 

Extending the infinitesimal scenario to finite, continu- 
ous processes in which the temporal evolution of p(t) and 
H(t) in time t is known, the mean heat and work can be 
found by integrating over the trajectory taken from p(0) 
and H(0) to p(r) and H(t), i.e. 

Q ■= [ dttv[p{t)H{t)] , (5) 
Jo 

W := I dttr[p{t)H(t)} , (6) 
Jo 

while the first law becomes 

AU:= f dt^-tr[p(t)H(t)} = Q + W . (7) 
Jo at 

The mathematical consistency of the above expressions 
and their compatibility with the predictions of thermo- 
dynamics have been verified for many models, for exam- 
ple, for processes that are induced by Markovian master 
equations [18]. 

There are two paradigmatic examples of all work and 
all heat processes that we introduce here and which will 
become important in the later part of the manuscript. 
The first process is a unitary process, which we will also 
refer to as closed, where the (non-equilibrium) evolution 
of the state is given by the Schrodinger equation, 

p{t) = - l -[H{t),p{t)]. (8) 
Mean heat and work are then 

Qumtary = ~ ~ / dt tl[[H (t) , p(t)} H (<)] = , 
W un itary = AU, (9) 



consistent with the physical intuition that no heat has 
been provided to the system during the evolution. The 
second example is a system that evolves through the ac- 
tion of a dissipative, i.e. open, Markov process via a 
master equation [26] , 

p(t) = -j i [H,p(t)}+C(p(t)), (10) 

with C being the dissipative Lindblad term. Under the 
assumption that the typical time scales associated with 
the time-independent H are much shorter than those as- 
sociated with C we can treat the system as almost iso- 
lated and use Eq.|l]) to compute its internal energy. In 
this limit Eq. ^ is valid with the Hamiltonian just being 
the time-independent H, 

dissipative 0, 

Qdissipative = tr[p(r) H] - tr[p(0) H] = AU, (11) 

which is in full agreement with the physical intuition that 
no work has been performed on the system. 

While these examples constitute special cases of con- 
tinuous processes the heat and work in a general process 
depend intimately on the exact details of the process. 
However, the caveat with this viewpoint is that in most 
real life applications one does not know what the dynam- 
ics of the state of the system is nor what the appropriate 
local Hamiltonian is at all times. Importantly, this is 
not just due to our ignorance of what happens at the 
quantum level. Quantum physics has strong fundamen- 
tal limitations on what we can know without choosing a 
measurement apparatus, measuring the system and inter- 
preting the data. Moreover, if the system is indeed mea- 
sured then the experimenter's choice of what degrees of 
freedom she actually measures will effect what the mea- 
sured heat and work will be. In other words, we propose 
that there is no one average heat and work for a partic- 
ular process, there are different sensible outcomes to this 
question and the answer depends on the choice of system 
Hamiltonians in time, H(t), that corresponds to specific 
measurement choices. 

The aim of this paper is to show that it is possible 
to formulate a general second law independently of these 
choices. To achieve this we will depart entirely from the 
traditional continuous trajectory approach and propose 
a rather drastic but pragmatic change of perspective. We 
develop a consistent framework of mean heat and work 
for discrete thermodynamic processes. The rationale for 
this approach is that while the true process is continuous, 
observations we make on the system are almost always 
discrete. (We will neglect here the possibility of moni- 
toring through continuous weak measurements.) For dis- 
crete snapshots of the dynamics, we find that by decom- 
posing the transition into possible sequences of two fun- 
damental primitives, it is possible to define heat and work 
for the discrete process in a way that is experimentally 
and mathematically clear. This allows us to establish a 
general second law for discrete processes between equi- 
librium and non-equilibrium states and the analysis of a 
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discrete Carnot cycle, where we uncover the usual Carnot 
efficiency. 

The paper is structured as follows. In Sec. [TT] we re- 
view the traditional perspective on the second law and 
the definition of entropy. In Sec. Ill we define the dy- 



namical configuration space of a system that allows us to 
formulate a notion of two primitives for discrete processes 
in Sec. |IV[ the discrete unitary and discrete thermalising 
transformations (DUTs and DTTs). Sec. [V] contains the 
main results of the paper. First we show that entropic 
inequalities when applied to discrete trajectories formed 
by concatenating DUTs and DTTs yields the second law 
of thermodynamics in the Clausius formulation. We then 
derive two consequences: We find the minimum and max- 
imum heat for a single DUT and DTT sequence and prove 
the existence of a discrete trajectory, formed by sequences 
of DUTs and DTTs, that connects two given thermody- 
namical configurations while asymptotically saturating 
the Clausius inequality. Finally we identify a discrete 
trajectory that connects the same initial and final config- 
urations as the continuous trajectory through a sequence 
of DUTs and DTTs, and which approximates the con- 
tinuous heat. In Sec. |VI| we derive the thermal efficiency 
of a discrete cycle, the Carnot efficiency, and conclude in 
Sec. Hill 



II. ENTROPY AND THE SECOND LAW 

In 1865 Clausius established that the overall heat flow 
in any cyclic, reversible process vanishes, implying that 
the integral over any non-cyclic process must be path 
independent. This led him to define the state function 
entropy, S, and the entropy change, AS, between the 
final and initial point of a reversible process, 



SQ 
T 



= 



T 



=: AS. 



(12) 



Clausius also showed that any cyclic process, reversible 
or irreversible, obeys 



SQ 
T 



< 0. 



(13) 



This relation is the basis for a formulation of the sec- 
ond law of thermodynamics, known as the Clausius- 
inequality. It is a statement for all thermodynamic pro- 
cesses, not just cyclic ones, that start from equilibrium 
at temperature T, 



< AS, 



(14) 



and it simplifies to Q < TAS when the system interacts 
with a bath at constant temperature, T. In this form 
Clausius' inequality establishes the existence of an upper 
bound to the heat received by the system. 

Clearly, Clausius's goal was to characterise different 
forms of energy and their interconversion. However, by 



formulating the second law of thermodynamics he de- 
fined a new quantity: entropy. In contrast, in modern 
information theory the focus is on the state of a system. 
Entropy is here used as the central physical quantity to 
measure the amount of information of a state, while heat 
and work, and energy in general, have no well-defined 
purpose for the interpretation of information processing. 
This opens the possibility of turning Clausius' original 
argument around! It allows us to use the entropy change 
in discrete quantum processes to define the average heat, 
and work. Before we proceed, let us first highlight that 
non-trivial entropy bounds exists for any process between 
two states. 

A state p describes an amount of information, quan- 
tified by the von Neumann entropy, S(p) = —tr\p hip]. 
The evolution of a quantum system from an initial state 
Pi to a final state pf through an arbitrary process, or 
quantum channel, has a meaningful associated entropy 
change, 



-tT[pflapf] +tr[/9 i lnp i ], (15) 



which quantifies the change of the encoded amount of in- 
formation. The entropy change is non-trivially bounded 
from above and below by virtue of the positivity of the 
relative entropy (classically Kullback-Leibler divergence 
[27]). The relative entropy, S(pi\\p2), between two states 
Pi and P2 characterises the number of additional bits re- 
quired to encode p\ when using the diagonal basis of p2, 
rather than the diagonal basis of p\. It is defined as [23] 

S(pi\\p2) ■= tr[pi lnpi] - tr[pi lnp 2 ] , (16) 

and is a positive quantity 

S(pi\\p2) > with equality iff p\= Pi- (17) 

Intuitively, the relative entropy is similar to a distance 
measure, however, it is important to keep in mind that 
it is asymmetric S(px\\p2) ^ S(p2\\pi)- Rewriting the 
entropy change, Eq. (15), in two ways 



AS{ Pl , Pf ) = -ti[Ap\np } ]+S{ Pl \\p f ) (18) 
= -tr[Ap In Pi ] -S{p f \\pi), (19) 

a lower and upper bound on the entropy change emerge 

- tr[Ap la > AS(pi,p f ) > -tr[Ap lnp f ]. (20) 

From the information theory point of view, bounds on 
the entropy change are important in their own right as 
they characterise how much information is lost or gained. 

If we now assume the special case that pf is a thermal 
state for the Hamiltonian Hf at an inverse temperature 
(3f then the lower bound becomes 



AS{pi,pf)>pftt[ApH f ]. 



(21) 



Interpreting tr[Ap7Jj] as the heat of the discrete process, 
the above expression would constitute exactly the second 
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FIG. 1: As a visual aid, points in dynamical configuration 
space are depicted in a (p, H )-coordinate system. A non- 
thermal configuration a = (pi, Hi) is shown as a blue square 
and a thermal configuration Cf(fif) = (pf,Hf)p, is shown 
as a red circle. The discrete trajectory from a to c/(/3/) is 
indicated by the black arrow. 



law of thermodynamics. This is exactly what we will 



pursue in Sec. IV e.g. in Eq. (28 1 



Interestingly, from ( 20 ) it is apparent that also an up- 



per bound on the entropy change exists that is rarely 
discussed in the literature. This maximum value of the 
entropy change is enforced to ensure that any reverse 
process, from pf to pi, also obeys the second law |32j . 



III. DYNAMICAL CONFIGURATION SPACE 

To assist our discussion of discrete quantum processes 
we introduce the concept of configuration space, follow- 
ing the spirit of [TBI HZ1 and propose a graphical 
representation for that space, see Fig. [T] 

Definition 1 Let S be the quantum system under inves- 
tigation, with Hs its Hilbert space, C(Hs) the set of lin- 
ear operators on Hs, and &(Hs) C C(Hs) the set of 
density matrices on Hs- We define the dynamical con- 
figuration space C(Hs) of S as the set formed by the pairs 
(p, H) = c with p £E 6(Hs) a density matrix \3ty and 
H G C(Hs) a Hermitian operator onHs whose spectrum 
is bounded from below. Points in the dynamical config- 
uration space c are called "configurations" to distinguish 
them from "states", p. 

The evolution of the system is described by discrete tra- 
jectories in C(Hs). 

Definition 2 A discrete trajectory T is defined as an or- 
dered list of elements of C(Hs) that describes the succes- 
sion of configurations, with each element (p, H) contain- 
ing both the density matrix p of S and the local Hamilto- 
nian H of S at that specific instance of the evolution. 

We stress that both p and H of a configuration point 
c € C(Hs) have a clear experimental meaning, p is the 
density matrix that one would reconstruct by state to- 
mography, i.e. the preparation of many copies of the 
same state p and the full tomographic measurement of 



its properties. The Hamiltonian of the system, H , is de- 
termined by the set of projective measurements {My}j 
the experimenter performs on the system to "measure 
the energy" together with the interpretation of the corre- 
sponding energy eigenvalues, Ej, so that H — J2j Ej Mj. 
(The choice of the measurement and interpretation can 
be motivated by a process tomography on the Hamilto- 
nian at any point in time. For this the system needs to 
be decoupled from the rest of the universe at that in- 
stance and evolve for a complete set of states for a short 
time interval r through the action of H. By measuring 
the final states of the evolution the unitary e~ lHr and 
hence H can be uncovered.) It is then straightforward 
to establish the internal energy and the entropy for each 
point in dynamical configuration space. 

Definition 3 For each configuration c = (p, H) in 
C(Hs) we define the internal energy as 

U(c)=tt[pH], (22) 

and the entropy as the von Neumann entropy S(p) of p, 



S(p) = -tr[p lnp]. 



(23) 



A central notion in thermodynamics is the canoni- 
cal Gibbs state, often also referred to as thermal state 
or equilibrium state. Since a thermal state, p, at tem- 
perature T is well-defined only with respect to a cer- 
tain Hamiltonian, H, it is actually the configuration 
c = (p, H) that is thermal. 

Definition 4 An element (p,H) E C(Hs) describes a 
thermal equilibrium configuration (or briefly thermal 
configuration^ if p is a Gibbs state of the Hamiltonian, 
H , for some finite inverse temperature (3 > 0, i.e. 



Z(P)' 



(24) 



with Z((3) = tr{e~@ H ] being the associated partition func- 
tion 13$ . In the following the thermal configurations will 
be indicated by c(/3) := (p,H)p with the subscript j3 spec- 
ifying the configuration's temperature. 

Thermal configurations (p,H)p are very special. 
Firstly, for a given Hamiltonian, H, from all possible 
states that have a fixed value of the internal mean energy, 
U = ti[pH], the thermal state maximises the entropy 
S(p) = — tr[plnp]. In other words c((3) — (p,H)p is the 
most unbiased configuration one can assert to the system 
given only the knowledge of U [30 . Another insight- 
ful characterisation of thermal configurations in terms 
of a property called complete passivity was achieved by 
Lenard [17] . building on ideas of Pusz and Woronow- 
icz |16j . Complete passivity captures the intuitive no- 
tion of thermal equilibrium. A configuration (p, H) is 
said to be passive if no work can be extracted from 
the system, i.e. W > cf. Eq. p}]), when subjected 
to any unitary transformation for a time r generated 
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by an arbitrary time-dependent Hamiltonian with the 
sole constraint that H(t) = H(0) = H. A configu- 
ration c = (p, H) is completely passive if all its regu- 
larised configurations c^ := {p® n ,H^ = £? =1 H j) 
are passive for n — 1,2,.... Here the unitary opera- 
tions entering in this definition are generated by arbi- 
trary time-dependent Hamiltonians H^ n \t) that satisfy 
the constraint H^(0) = ff(")(r) = H (n) . I.e. during 
t s]0, t[, (i) is allowed to introduce any sort of inter- 
actions between the various copies of p. It turns out that 
while all c = (p, H) with commuting p and H are passive 
configurations, only thermal configurations c(/3) and the 
ground state are completely passive [THlIIZj- To stress the 
special role of thermal configurations graphically, they 
are denoted as red circles while all configuration that are 
not thermal will be called non- equilibrium configurations 
and are depicted as blue squares, see Fig. [TJ 

A note on gauge. Given a generic state p € @(Hs) 
which is full rank, there always exists a Hermitian opera- 
tor H € C(7is) an d a P > such that (p, H)p is thermal 
[55] at the inverse temperature j3. In fact the problem 
admits infinite solutions, since there are two gauge free- 
doms for the choice of H and f3. Firstly, the zero-point of 
the energy scale can be chosen arbitrarily by a constant 
a. The second gauge, 6, is the temperature itself, which 
sets a spacing of the energy scale. The pair {H, 0} is 
equivalent to {b(H + a),/3/b} in that they have the same 
set of thermal configurations. In particular the internal 
energy ( 22 ) and the entropy ( 23 ) of such configurations 
do not depend on the values of a and b. In the follow- 
ing we will assume that both gauges have been chosen to 
some fixed values. 



IV. DISCRETE TRANSFORMATIONS IN 
DYNAMICAL CONFIGURATION SPACE 

Among all possible discrete transformations in dynam- 
ical configuration space C('Hs) we identify two classes 
that admit a clear analysis of the energetic balance and 
can be used as primitives for general discrete dynamical 
evolutions. 



A. Discrete Unitary Transformations (DUTs) 

These transformations map an initial configuration 
c i = (Pi)Hi) to a final configuration Cf = (pf,Ht), de- 

with the only constraint that 



noted as Cj ^> c / 



H 




Pi 



P2 P3 



FIG. 2: Originating from the initial state Ci — (pi,Hi) three 
DUTs are shown, depicted as blue arrows, each ending at 
a final configuration Cj = (pj = Vj pi Vj , Hz) for j = 1, 2, 3. 
The special case where Cj is transformed into c 1 is an example 
of a DUQ. Also shown are two DTTs, depicted as red fuzzy 
arrows, originating from a and ending at thermal configura- 
tions ci(/3i) = (p 1 ,H 1 ) l3l and c 2 (^2) = (p2,Hi)/s 3 where /3i 
and fa are inverse temperatures associated with the thermal 
configurations. 



arbitrary unitary evolutions V — Te 1 dtH (t)/ h where 
T indicates time ordering. For the discrete mapping 

Ci cf no assumption is made on the time duration r 
nor the specific form of H(t) which realises the unitary 
V. In analogy to the continuous situation, we define the 
work done on the system due to a DUT identical to the 
total variation of the internal energy, At/, i.e. 



W( Cl D ^ T c f ):=U(c f )-U(c t ) 

= tv\pi (V^HfV-Hi 

while no heat is associated with DUTs, i.e. 



dut . 
Q(Ci -> Cf) 







(26) 



(27) 



Pf = V Pi ^ 



(25) 



A special class of DUTs are the Discrete Unitary 
Quenches (DUQs). Experimentally, a quench is an 
abrupt, instantaneous change of the system Hamiltonian 
which leaves the system density matrix unchanged, i.e. 

V = 1, (pi,Hj) (pi,Hf). We also note that for full 
rank states pi a DUQ can be found that brings (p i; Hi) to 
a final configuration that is thermal, c/(/3) = (pi,Hf)p, 
with the Hamiltonian defined as Hf — — i(lnpj + InZ). 

DUTs will be denoted as blue arrows in the graphi- 
cal representation of the configuration space, see Fig. [2] 
These "work arrows" indicate the closed nature of the 
transformation. The DUTs characteristic properties are 
summarised here: 



for some unitary V. No constraint is posed on the rela- 
tionship between Hi and Hf. The definition of DUTs 
is inspired by continuous unitary transformations, see 
Eq. ([8]). There the system is thermally isolated while 
evolving through the action of some external force that 
modifies the Hamiltonian in time, H(t), and generates 



1. DUTs can start from and end at either non-thermal 
or thermal configurations. 

2. DUTs change the state by a unitary and no entropy 

DUT 

change is induced by a DUT, i.e. AS'(c i — > c/) = 
S(pf) - SQh) = 0. 
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3. DUTs can be concatenated to produce another 
DUT. 

4. Any DUT has an inverse that is also a DUT. 



B. Discrete Thermalising Transformations (DTTs) 

DTTs are defined as those transformations which take 
a generic Cj = (pi , H) into a Gibbs state at some in- 
verse temperature j3, c/(/3) = (pi = e z , H^j , without 

modifying the system Hamiltonian H. The prototypical 
example of a DTT is an arbitrary thermalisation process 
in which the system is put into a weak thermal contact 
with a reservoir at inverse temperature /3 and left un- 
til its state becomes time- independent. Physically this is 
realised by the system weakly interacting with a large ex- 
ternal environment. The requirement of a small coupling 
ensures a clear definition of a local system Hamiltonian. 
For example, the dissipative evolution (p(t),H) defined 
in Eq. ( 10 ) with the additional assumption that the Lind- 



blad term C commutes with H will for t — > oo converge 
to pi. In analogy with this continuous process, we as- 
sume that the internal energy change due to a DTT is a 
result solely of the heat absorbed by the system 



Wc f (l3)):=U(c f (P))-U( Cl ) 



= tT[H( Pf - Pi )], 

while the work of a DTT vanishes, 

W(c,% T cf{fi)) :=0. 



(28) 



(29) 



This non-trivial expression of the heat is exactly of the 
form that we expected from the bounds on the entropy 
in Eq. (21 ), and it will be the basis for deriving a general 



second law for discrete quantum trajectories in the next 
section. 

DTTs will be denoted as horizontal red arrows in the 
graphical representation, see Fig. [2| The fuzziness of 
these "heat arrows" indicates the open nature of the 
transformation. The characteristic properties of DTTs 
are summarised here: 

1. DTTs always end in thermal configurations. 

2. DTTs do not change the Hamiltonian. 

3. The entropy change associated with a DTT is in 

DTT 

general nonzero, i.e. A5(c, — > c/(/3)) = S(pf) — 
S( Pi ) + 0. 

4. DTTs can be concatenated to produce another 
DTT. 

5. The inverse of a DTTs is in general not a DTT. 
Only if the initial configuration Ci was already ther- 
mal, can the action of a DTT be reversed by an- 
other DTT. 



V. HEAT AND CLAUSIUS INEQUALITY FOR 
DISCRETE THERMODYNAMIC PROCESSES 

Having identified two fundamental process primitives 
in configuration space, we now focus on more complex 
discrete trajectories. These can start from equilibrium 
or non-equilibrium configurations, however, we restrict 
ourselves to discrete trajectories that can be obtained by 
concatenating DUT and DTTs. Within this scenario we 
will be able to formulate a general second law for dis- 
crete quantum processes, that does not require detailed 
knowledge of the continuous state and local Hamiltonian 
evolution. 



A. Single DUT+DTT transformations 

Let us begin with the simplest non trivial discrete 
transformation which can be used to connect two equi- 
librium configurations. 



1. Equilibrium to equilibrium processes 

We consider a trajectory that starts from a thermal 
configuration Cj(ft) = (pi 1 Hi)p i and ends at a final ther- 
mal configuration c/(/3/) = (/?/, Hf)p, via the action of a 
single DUT followed by a DTT. The heat of the discrete 
process can then be determined as the sum of the heats 
of each component, for both of which the heat is a well- 
defined quantity. The DUT first unitarily rotates the 
input density matrix to pi — V pi while the Hamilto- 
nian changes from Hi to Hf, ending in an intermediate 
(not necessarily thermal) configuration c\ — (pi,Hf). A 
DTT follows that brings ci to e/(/3/), resulting in the 
discrete overall trajectory 



DUT DTT , n , 
— > c l — > c /(P/)> 



(30) 



shown in Fig. [3j While heat of process ( 30 ) is only ex- 
changed during the DTT, the amount of exchanged heat 
depends on the DUTs unitary V 



Q( Cl u ^ c f {pf)) = tr[( P/ -Vpi V*) H f ] 



(31) 



Clearly, the value of the heat depends on the choice of the 
unitary V with the maximum and minimum heat given 
by 



Qt 



N 

E^/( fc )( ! 



fe=i 

JV 



-0fH f (k) 
-PiHf(k) 



-PiHt{k) 



-PiHi(N-k+\) 



k=l 



Z 



f 



(32) 



where {Hf(k)}k and {Hi(k)}k are the eigenvalues of Hf 
and Hi ordered in decreasing order and Z it f the parti- 
tion functions of the initial and final configuration. The 
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H 
H f 

Hi 



Cl 



■-^• c/ (^) 




Pi 



pf 



FIG. 3: Trajectory connecting two Gibbs configurations, 
Ci{Pi) = (Pi,Hi)fj z -t Cf(f)f) = (pf,Bf)pt (black arrow) and 
discrete decomposition into a DUT, Ci(pi) — > ci (blue arrow), 
to an intermediate point ci = (pi,Hf) followed by a DTT, 
ci — > Cf(pf) (red arrow). 



derivation of this expression and the corresponding uni- 
taries V m i n and Umax, are given in Appendix [A"| 

However, for any possible choice of the DUT connect- 
ing Ci(j3i) to the intermediate step ci, i.e. for any uni- 
tary transformation V, a second law can be established 
by linking the heat and the entropy change. The entropy 
change is bounded according to Eq. ( 20 1 , by 



H 
Hf 

H„ 
H, 

b 

H 

H { 



■ — > • 

■ — -^t 

^< c m (An ) 

» 

cm 



■-»■• 

t 



Cl(ft 



&S(pi,Pf) = S(pf) - S(pi) 

> -tr[(p/ - pi) hip/] 

= /3 / tr[(p / -7 P iW) J ff / ]. 



implying 



(33) 



(34) 



Thus the process ( 30 1 obeys a Clausius-type inequality, 



Pf 



FIG. 4: Panel a shows the realisation of a general trans- 
formation Ci(/3i) — > Cf(f3f) as a sequence of two DUT+DTT 
transformations, via an intermediate thermal configuration, 
CmiPm) = {pm,H m )p m . Panel b shows the realisation of a 
general transformation Ci{fii) — > c /(/3/) through a sequence 
of DUT+DTT transformations, via many intermediate ther- 
mal configurations, c k (/3 k ) = (p h ,H k )p h . 



allows us to identify the heat associated with T with the 
quantity 

Q(T) = Q(ci °^ c 2 (/? 2 )) = tT[( Pf -V Pi V*)Hi]. (36) 



cf. ([14], which states that the heat absorbed by the sys- T h e lower bound of Eq. (f20b can then be used to estab 



tern is upper bounded by the entropy change. 



2. N on- equilibrium to non- equilibrium processes 

We now turn to discrete non-equilibrium processes for 
which establishing the Clausius inequality in the contin- 
uous case has only recently been addressed [35]. In our 
approach this can be done by observing that given two 
generic configurations c; = (pi,Hi) and c/ = (pf,Hf) in 
C(Hs), it always possible to connect them via a discrete 
trajectory composed by three primitive steps which dif- 
fers from the one given in Eq. ( 30 1 only by a final DUQ 
transformation. Specifically we can write 



lish a Clausius-type inequality for the discrete transfor- 
mation (|35|), i.e. 



AS( Pl ,p f ) = AS( Pl ,p 2 ) (37) 
> -tr[(p f -Vp i V^bxp f ]=0 2 Q(r), 

where /3 2 is the temperature of the intermediate configu- 
ration C 2 (/? 2 ). 



B. Sequences of DUT+DTTs 



The trajectories defined in Eq. (30) and Eq. (35) 



DVT DTT to \ DUQ 

T := c< — ► Ci — > C2{p 2 ) — > c f , 



(35) 



with intermediate configurations c\ = (V piV',H\) and 
c 2 (/3 2 ) = [pf-,Hi)fe- Notice that the first and the last 
step of ( [35] ) do not alter the entropy of the system, nor 
contribute to the heat exchange since they are DUTs. 
This implies the identity AS(pi, pf) — AS l (pi,p 2 ) and 



are just specific choices of discrete trajectories connect- 
ing two configurations Cj and c/. We will now show 
that a Clausius inequality, e.g. inequalities of the type 
(34), holds for general discrete processes as long as 



they can be decomposed in a sequence of DUT+DTTs 
steps. To show this, we first consider the discrete tra- 
jectory, 7, pictured in panel a of Fig. H where Cj(/3j) = 
(pi,Hi)^ i is transformed into c/(/3/) = (pf, Hf)p, via two 
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DUT+DTT transformations and a third thermal config- 
uration c m (f3 m ) = (Pm,H m )j3 m , i.e. 

,„s DUT+DTT (a v DUT+DTT (Q s , QtA 

7 := Ci(/3i) — > Cm{p m ) — > CfiPf). (38) 

In this scenario the following inequality for the entropy 
holds, 



AS(pi,Pf) = AS(pi,p m ) + AS(p m ,Pf) 

^ Pm Qm + Pf Qf, 



(39) 



is obtained by summing over the A contributions of the 
various steps, i.e. 



N-l 



A(T) := J] Mck(Pk) 



DUT+DTT 



Cfc+x^fc+i)), (42) 



fc=0 



with fc = and fc = N corresponding to the initial and 
final configurations, i and /, respectively. Then the lower 
entropy bound, Eq. (40), can be neatly expressed: 



DUT+DTT 



where Q m = Q(c;(/3j) 

rM la \ DUT+DTT (a ^ 

used for the two DUT+DTT transformations 

Inequality p9l can immediately be generalised to 



c m {f3 m )) and Q f = 
and where Eq. (34) was 



an arbitrary number of intermediate DUT+DTT steps 
connecting Ci(/3i) to c/(/3/). Specifically, consider a 
generic discrete trajectory T composed of N consecutive 
DUT+DTT steps that pass through the thermal config- 
urations {Ci(Pi), Ci(Pi), 02(02), •" ,CN-l(0N-l)>Cf(Pf)} 

as shown in panel b in Fig. [4] Then by expressing the 
total entropy increment AS(pi 1 pj) as a sum of terms 
AS(pk, Pk+i) associated with the various steps of T and 
applying Eq. (34) to each one of them, the Clausius in- 



equality becomes 



AT-l 



AS( Pi ,p f )> Y,p- k +iQ(c k (p k ) 



DUT+DTT 



fe=0 



Cfe+l(^fe+l))- 

(40) 



The generality of this derivation implies that sequences 
of discrete unitary and discrete thermalising transforma- 
tions always fulfil a Clausius type equation. 

To formulate this as a lemma, we introduce a useful 
discrete process quantity, A, for a DUT+DTT sequence, 



. f f a s DUT+DTT f a n 

H c k(Pk) — > c k +i(Pk+i)) 

V k p k vl) H k+1 ], 



'■— Pk + l^[(fik + l - •'fcP/c ^fc 

a n( <a \ dut+dtt 



(41) 



c k +i{P k +i)). 



The quantity A is the discrete analog to the continuous 
expression J For the trajectory T the overall A (7~) 



Lemma 1 Any trajectory T made of sequences of 
DUT+DTTs fulfils a Clausius inequality of the form 

AS( Pi ,p f )>k(T). (43) 
While in general a single DUT+DTT process cannot 
saturate the equality, see Appendix |XJ we will now show 
that augmenting intermediate steps will always increase 
the r.h.s. of Eq. (40). Moreover, we find that in the limit 



of infinitely long sequences the asymptotic saturation of 
the inequality (40) is always possible. 



C. Saturating the Clausius bound 



To show that the entropy bound Eq. (40) can be sat 



urated we construct a class of trajectories T' from a 
generic DUT+DTT trajectory T, as depicted in panel 
b of Fig |4j for which the functional A(T') is always 
larger than A(T). One class of trajectories T' is the 

DTT 

trajectory identical to T however with the step c k — > 

DUQ+DTT 



Ck 



Cfc+i(At+i) replaced with the sequence 
c m (Pm) DU 9±f TT Ck+i(Pk+i) as shown in Fig. [5J where 
Cm(Pm) = (Pm,H m )p m has an intermediate density ma- 
trix 



Prn =PPk + (l~ P) Pk+1 



(44) 



with mixing probability p g]0, 1[ |36) . The increment on 
the r.h.s. of (43) for the new trajectory T', A(T',T) = 
A(r')-A(T) = ^ OT Qi+ ( a fe +iQ 2 -/3fe+iQfe, see Fig.g is 
then strictly positive for any p 



A(T',T) = -tr[lnp m (p m - p k )} - tr[lnp fe+ i (p k +x - Pm)} + tr[lnp k+1 (p k +i - p k )\ 

> -pS(pk\\Pk) - (! ~P) S(p k +i\\pk) - pS(p k +i\\pk) -(l-p) S(p k +i\\pk+i) + S(p k+ i\\pk) 

-0, 



where we have assumed p k 7^ Pk+x, an d used Eq. (44) 
and the joint convexity of the relative entropy [2"B] . 



with equality iff p k +i = pi = p k . We summarise this 
result in the following Lemma: 



S(pi\\pPk + (1 —p)pk+i) < pS(pi\\p k ) + (1 - P)S(pi\\pk+i) Lemma 2 Adding intermediate thermal configurations 
S(pp k + (1 -p)pk+x\\pi) < pS(p k \\pi) + (1 - p)S(p k+l \\pi) 

(45) 
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H 




Ck+l(Pk+l) 



Pm 



Pk+l 



FIG. 5: Discrete trajectory connecting Ck — > Cfc+x(/3fe+l) 
by a single DTT with heat Qfc, and by a sequence of two 
DUQ+DTT transformations via an intermediate thermal con- 
figurations, c m (/3 m ) = (p m , H m )/3 m . The two DTTs have 
heats Qi and Q2, respectively. 



c m (/?m) (see Eq. 44) to an U trajectory T results in a new 



trajectory T' with increased A, i.e. 

A(T',T) = A(T') - A(T) > 0. 



(46) 



Having confirmed that it is possible for any given dis- 
crete trajectory to introduce intermediate steps such that 
the r.h.s. of Eq. (43 1 increases, the task is now to show 



that the entropy bounds can be saturated by reiterating 
the procedure. The proof relies on lower bounding A and 
showing that the bound converges to the upper bound 
on A, Eq. (43), in the limit of infinite steps. The detailed 



derivation is given in Appendix [B] proving the following 
theorem: 

Theorem 1 Let T be a discrete trajectory connecting the 
initial Gibbs configuration Ci{f3i) = (pi,Hi)p i to the final 
Gibbs configuration Cf(f3f) = (pf,Hf)p, via a sequence 

of N concatenated DUT+DTT steps Ck{fik) Dur £^-? TT 
Cfc+i (/3fc+ 1) connecting the thermal configurations T = 
{cj(ft) = c (/3o) ) ci( ( 8i),...czv(/3iv) = Cf(/3 f )} as in panel 
b in Fig. [^J Then a sequence of trajectories 7^ exists, 
obtained from T by adding n intermediate thermal steps, 
which saturates the Clausius bound (40) in the asymptotic 
limit, i.e. 



AS( Pi , Pf )= lim A(X). 



(47) 



D. 



Approximation of continuous processes by 
discrete processes 



In the introduction we have seen that for continu- 
ous processes where consistent definitions of p(t) and 
of the local Hamiltonian H(t) can be assigned for all 
t, Eq. ([5]) defines the heat absorbed by the system. We 
have already discussed the difficulties of knowing p(t) and 
identifying a proper local Hamiltonian H(t) for the sys- 
tem. However, in what follows we will assume that some 
"valid" continuous trajectory c(t) = (p(t),H(t)) G C{% s ) 



is given for which Eqs. ^ and (|6| apply. We now wish 
to identify a discrete trajectory that connects the same 
initial and final configurations as the continuous trajec- 
tory through a sequence of DUTs and DTTs, and which 
approximates the continuous heat. The analysis leads to 
the following theorem: 

Theorem 2 For a continuous process between two con- 
figurations Co and c T that obeys the Clausius inequality, a 
discrete trajectory V exists that connects the same config- 
urations and has exactly the same heat as the continuous 
process. 

Proof: Consider an infinitesimal heat increment along 
the continuous trajectory, 



6Q(t)=tr[(p(t)-p(t-dt))H(t)}, 



(48) 



with p(t) and p(t — dt) being the density matrices of two 
infinitesimally separated configurations on the trajectory. 
Define the initial and final configuration for a discrete 
trajectory to be C, = (p(t — dt),H(t — dt)) =: (p±,Hi) 
and cj = (p(t),H(t)) := (p 2 ,H 2 ). To compare the con- 
tinuous heat ( [48] ) with a discrete heat we need to identify 
a discrete trajectory, 7, connecting the same initial and 
final configuration as the continuous trajectory. 

One example is the sequence 7 shown as a solid line in 
Fig.© 
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DUQ DTT - DUQ „ DTT / -t \ DUQ 
Ci -> Ci -> c 2 (l) c 2 c 2 (l) C/, 



(49) 

where c 2 (l) = (/5 2 , H 2 )p=i and c 2 (l) = (p2,H 2 )p=i are 
equilibrium configurations |37j and c 2 = (/5 2 ,-ff 2 ) is a 
non-equilibrium configuration. The inverse temperatures 
of the equilibrium configurations c 2 (l) and c 2 (l) are both 
chosen /3 = 1. Heat Qi, 



, DUQ+DTT _ r, ( DTT-,,, 

Qi: = Q(Ci — > c 2 (l)) = Q(ci — > c 2 (l)J 
= -tr[(p 2 - pi)lnp 2 ], 



(50) 



is exchanged when passing from cq — > c 2 (l) via a 
DUQ+DTT through the intermediate configuration ci, 
see Fig. [6] Heat Q 2 , 

DTT 

Qi ■■= Q(c 2 c 2 (l)) = -tr[(p 2 - p 2 ) In pa], (51) 
is exchanged when passing from c 2 — > c 2 (l). Therefore 

Q(rr) = Q t + Qi- (52) 

On the other hand, the continuous heat increment 
SQ(t) can be decomposed into two heat contributions, 

SQ(t) = tr[(p(*) - p(t - dt))H(t)] = tr[(p 2 - Pl )H 2 ] 

= tr[(p 2 - P2)H 2 ] + tr[(p 2 - pi)H 2 ] 

= -tr[(p 2 - p 2 )lnp 2 ] - tr[(p 2 - pi)lnp 2 ] 

= -Qf + Qi- 

(53) 
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FIG. 6: The figure shows the initial and final non-equilibrium 
configurations, c,: and Cf, of the continuous process. A dis- 
crete trajectory, 7, that connects Ci and c/ is also indicated 
(solid). Also shown is the closed loop discrete trajectory u> 
(dotted). The two thermal configurations £2(1) = {f>2,H2)p=\ 
and 02(1) = (p2, H2)p=i are both chosen with inverse temper- 
ature /3 = 1, as indicated in the circles. 



where Q / would be the heat absorbed by the system if it 
passed from c/ to £2(1) via a DTT, i.e. 

Qf ■= Q(c f c 2 (l)) = -tr[(p 2 - pa) lnp 2 ]. (54) 
To compare the continuous heat |53|) with the discrete 



FIG. 7: Representation of the cycle of Eq. ( 59 I 



heat (52), we introduce a second discrete trajectory ui. 



This is a closed loop sequence of DUQ and DTT trans- 
formations, see Fig. [6j 



„ MS DUQ* DTT , , DU Q DTT-,., 

cj:=c 2 (1) -> c 2 -> c 2 (l) -> c/ -> c 2 (l). 



In trajectory ui heat is exchanged from c 2 — ^ c 2 (l) and 
from Cf —J- c 2 (l), so Q{uj) = Q2 + Qf- As discussed 
in previous sections the discrete heat always obeys the 
Clausius inequality ([34]), and with /3 = 1 for both steps 
this implies 



= AS(p 2 ,p 2 ) >Q 2 + Q f 



(56) 



for trajectory uj. Using this in Eq. (52) we find that the 
heat associated with the discrete trajectory 7 from c$ to 
Cf is a lower bound to the infinitesimal continuous heat 



for the same initial and final configuration, (48), i.c 



Q(rf) = Qi + Qi <Qi-Q f = SQ(t) 



(57) 



This can immediately be extended to the full continuous 
process: For any arbitrary continuous process between cq 
and c T there is always a discrete trajectory T between the 
same two configurations that has a lower heat than the 
continuous heat, Eq. ([5]). Moreover, by augmenting inter- 
mediate steps in the discrete trajectory T, resulting in the 
trajectory V that passes through an infinite sequence of 
points c(t), it is possible to increase the associated heat, 
as shown in Sec. lVBl From Theorem [T] follows that if the 
continuous trajectory fulfils the Clausius inequality, then 
a discrete trajectory connecting the same initial and final 
configuration can be found that has the same heat as the 
continuous trajectory. ■ 



VI. THERMAL EFFICIENCY 



The last piece in our analysis of the energy balance in 
discrete quantum processes is to determine the efficiency 
of a discrete cyclic process, such as the one depicted in 
Fig. 7|where c x (^i) = {px,H x )p x , c 2 (/3 2 ) = {p 2 ,H 2 )p 2 are 
equilibrium configurations, while C3 = (pi,H 2 ), and C4 = 
(p 2l Hi) are not. This results in the following Lemma: 



Lemma 3 The thermal efficiency of the discrete cycle 
depicted in Fig. [?] is bounded by the Carnot efficiency, 



(58) 



(55) and the optimal efficiency is achievable. 



Proof: For the overall loop, 



,„ \ DVT DTT .a \ DUT DTT to \ /cr^ 

ci(p\) -> c 3 -> c 2 ((3 2 ) -> c 4 ci(ySi) , (59) 



the entropy change is zero. However, the entropy change 
nevertheless bounds the heat of the two heat-producing 
DTT processes, C3 — > c 2 {(i 2 ) and C4 — > ei(p\), 



AS = ^ h Q{c 3 -> c 2 ) + f3i Q(a -> ci). 



(60) 



This implies that at least one of the two heats must 
be negative. Let us assume for instance that the heat 
exchanged with the thermal reservoir at temperature 
T 2 = l/(fcs/3 2 ) is positive, Q(c s — > c 2 ) > 0, while the 
other heat is negative Q(c4 — > c\) < (other scenarios 
can be treated analogously, see below). The total heat 
absorbed per cycle is 



Q(ci -> ci) = Q(c 3 -> c 2 ) + Q(c 4 -> ci), 



(61) 



with the energy balance implying that the overall ab- 
sorbed heat must be equal to the negative work done on 
the system during the cycle, 

= AU(a ->■ ci) = Q( Cl ci) + W{ci ci). (62) 
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The thermal efficiency is defined as the ratio between the 
work performed and the heat absorbed, leading to 



-W(d -t ci) Q(ci -> ci) 



Q(c 3 -> c 2 ) Q(c 3 -> c 2 ) 
_ Q(c 3 -> c 2 ) + Q(c 4 -> Ci) 
Q(c 3 -> c 2 ) 

1 - — = 1 £ , 



(63) 



where we used Eq. (60). If T\ ^ T 2 the system absorbs 



heat from a higher temperature bath and gives heat to a 
lower temperature bath. The efficiency ry is then positive 
and smaller than unity with the optimal efficiency repro- 
ducing the Carnot efficiency. The optimal efficiency can 
be reached by augmenting the discrete trajectory to satu- 
rate the equality in the Clausius-inequality, see Theorem 

HI ■ 

Remark: If instead T 2 ^ 7\, i.e. heat is absorbed 
from a lower temperature bath and given to one at a 
higher temperature, the system operates as a refrigera- 
tor. In this case the total work absorbed by the system, 
W{c x -> ex) = Q(c 3 -> c 2 ){T 1 /T 2 - 1), is positive. The 
efficiency of the process can be measured by the coef- 
ficient of performance, COP coo /m g , defined as the ratio 
between the heat absorbed from the cold reservoir T 2 (i.e. 
Q(c 3 — > c 2 )) and the total work done on the system 



COP 



cooling 



Q(C3 



C2 



T 9 



(64) 



which again is always smaller than one. Finally, if the 
signs for the heats Q(c 3 — » c 2 ) and Q(c± — > c%) are in- 
terchanged the above argument still holds with Eqs. ( [63] ) 
and (64) being replaced by the inequalities r) < 1 — 

Ti 



and COP 



cooling ^ 



-Ti ' 



respectively. 



VII. CONCLUSIONS 

The early development of thermodynamics culminat- 
ing in the formulation of the second law also gave birth to 
a new quantity, the entropy, whose physical meaning was 
at first opaque. Only later was its meaning elucidated 
by the works of Boltzmann and others. In this paper we 
proposed to turn the original argument around and use 
the well-established notion of entropy that characterises 
the information content in a (quantum) state to motivate 
the definition of a notion of heat for discrete quantum 
processes. The approach circumvents the large cluster 
of problems surrounding the idea of a unique definition 
of heat and work in processes where (i) the Hamiltonian 
of the system is not well-defined due to the open nature 
of the system and (ii) there arc fundamental limitations 
on the notion of trajectories where full knowledge is only 
given at discrete points in time when a measurement with 
a specific Hamiltonian occurred. 



By introducing thermodynamic configurations, iden- 
tifying two primitives for discrete processes, DUTs and 
DTTs, and defining heat to pertain only to DTTs we 
were able to uncover a general second law valid for any 
discrete process consisting of sequences of DUT+DTTs 
between both, equilibrium and non-equilibrium configu- 
rations. Moreover, we showed that an infinite sequence 
of DUT+DTT processes exists that saturates the Clau- 
sius inequality. In other words, saturation occurs when a 
discrete trajectory is mapped out into a continuous one 
by a sequence of measurements that are infinitely close 
together. This provides a link between reversibility - here 
the reversibility of a discrete process - and the equality 
in the second law for discrete processes, reminiscent of 
the Clausius' statement of equality for thermodynamic 
reversible continuous processes, Eq. (12). On the other 



hand, we also showed that for any continuous process 
between two configurations that obeys the Clausius in- 
equality, there exists a discrete process between the same 
configurations with the same heat. Finally, we showed 
that for the discrete version of a thermodynamic cycle, 
formed by a discrete trajectory passing through four con- 
figurations Carnot's efficiency is recovered. 

The strength of our approach is to give meaning to heat 
and work, reversibility and efficiency following from just 
a few sensible and simple definitions. In some respect 
this is analogous to the axiomatic approach to thermo- 
dynamics first developed by Caratheodory [31 . We hope 
that the presented analysis will inspire discussions and 
future work on characterizing heat and work in quan- 
tum processes. Of course many open questions remain. 
One direction of particular relevance is clearly the iden- 
tification of a proper metric in configuration space, that 
would allow to quantify, in a precise and (hopefully) op- 
erationally well defined way, how distant two generic dis- 
crete trajectories are. 
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Appendix A: Minimum and maximum heat for a 
single DUT+DTT process 



Here we discuss the impact of the unitary V of Eq. ( 31 ) 
on the heat of the DUT+DTT process pob. Specifically, 



we want to identify the DUTs that maximize and min- 
imize the heat on the r.h.s of Eq. (34) and establish if 



a DUT exists that leads to a saturated Clausius equal- 
ity. The last question can be easily solved by noticing 
that saturation occurs only when the inequality holds in 
Eq. (33). However, this is only true iff p\ = pf, i.e. pi 



and pf must be unitarily equivalent for some unitary Vq, 
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FIG. 8: Realisations of the map from the Gibbs state Ci{j3i) = 
{pi,Hi)p i to the Gibbs state c/(/3/) = (pf,Hf)p. via some 
DUT, specified by a unitary V , followed by an appropriate 
DTT. 



VoPiVg = pf. In other words equivalence in (30) requires 
that no DTT enters in the process so that ASljH, pf) =0. 
For any non-trivial DTT a finite gap between the l.h.s. 
and the r.h.s. of Eq. (34) exists. (This is not true how- 
ever for sequences of DUT-DTT transformations as con- 



sidered in Sec. V C 



To determine the minimum/maximum gap we require 
the maximum/minimum heat with respect to all possi- 
ble DUTs of a DUT+DTT process connecting two ther- 
mal configurations Cj(/3j) = (pi,Hi)p. and Cf(/3f) = 
(Pf' H f)f3fi see Fi g[§ i-e. 



Qmax := maxtr[(p/ - V pi V r ) Hf], 



mmtr[(pf-V Pi V*)Hf}, 



(Al) 



where the maximisation/minimisation is taken over all 
unitary transformations V. This task is solved with the 
following Lemma. 

Lemma 4 Given A = J2j a j\ a j)( a j\ an d B = 
Pj\j3j){Pj\ Hermitian operators on a finite dimen- 
sional Hilbert space %, with otj and fjj being their cor- 
responding eigenvalues which are ordered in decreasing 
order (i.e. ctj ^ ctj+\, fij ? fij+i)- Then the minimum 
of tr[AV B V^} over the set of unitary transformations is 
achieved by the unitary V m i n which maps the eigenvector 
{\/3j)} of B into the eigenvectors {\ctj)} of A in such a 
way that 



KiinlA,) = \0tN-j + l)- 



(A2) 



i.e. the maximum eigenvector of B is mapped into the 
minimum eigenvector of A. As a consequence the mini- 
mum expectation value is 

min tr[AVBV*] = ti[A V min B V^J = ^ a^ N ^ +1 . 

(A3) 

Similarly the maximum of tv[AV BV^\ over the set of 
unitary transformations is achieved by the unitary V maiyi 



that maps the eigenvectors {|/3j)} of B into the eigenvec- 
tors {\oij)} of A in such a way that 



(A4) 



Consequently the maximum expectation value is 
maxtrL4 VBV*]= tr[A V max B V^J = ^ atfy 

(A5) 

Proof: These minimum and maximum expectation values 
are a trivial consequence of the Theorem 2 of Ref . [T7] . ■ 



Application of Eqs. (A3) and (A5) gives the minimum 



and maximum heat for the DUT 



-DTT process (30) 



N 

Qmax = ]Ttff(fc)( 

fc=l 
N 

Qmi„ = ^^/(fc)(^ : 

fc=l ^ 



Zf Zi 

-PfH f {k) p -PiHi{N-k+l) 



Zi 



(A6) 

where {Hf(k)} k and {Hi(k)} k are the eigenvalues of Hf 
and Hi ordered in decreasing order. 



Appendix B: Proof of Theorem \T\ 

Consider the fc-th step of the trajectory T in panel b 
in Fig. [4j which connects the thermal points Cfc(/3fc) = 
(p k ,H k )p k and c k+1 (/3 k+1 ) = (p k+ i, H k+1 )p k+1 . We de- 
fine a new trajectory Tl. n which is identical with the 
original trajectory T except that the fc-th step is now 
replaced with a sequence of n — 1 intermediate thermal 
configurations c[(l), • • • , c' n _ 1 (l) that are linked through 
a DUQ-DTT sequence 



DUQ+DTT , 



DUQ+DTT , DUQ+DTT 



Cfc+1- 
( B1 ) 

The configurations c' e (l) have density matrices 
p' n _l defined by the mixtures 



P't = 1 - 



Pk 



Pk+1 for I = 0, • • • , n, (B2) 



where p' = p k and p' n — pk+i- Applying th e C lausius 
inequality ( 40 1 in the form of ( 43 ) to sequence ( Bl ) yields 



AS( Pk , Pk+1 ) > k{Ti n ). 



(B3) 



By definition (41), the transformations being discrete 
unitaries and the fact that intermediate inverse temper- 
atures are all set to 1, A(7^. n ) can be expressed as 

n-l 



A(T fc '.J = X>[(/^ - p' e+1 ) \n P ' e+1 ] 



£=0 



= tr 



(Pk 



Pk+ 



^ n— 1 

1) -J2 ln ^+i 



£=0 



(B4) 
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On the other hand the entropy change AS(p k , p k +i) can 
be lower bounded by 



n-l 
1=0 

= tr 



n-l 



(p fe - p k+ i) - V lnp^ 



A(7;' )+tr 



fc - pk+i) 



lnp fc - lnpfc+i 



= A(7]L) + 



^(pfe+illpfe) + S(p k \\pk+i) 



which implies 
A(T fc '.„) > AS{ Pk , Pk+1 )- 



(B5) 

S(p k+1 \\p k ) + S(p k \\p k+1 ) 
n 

(B6) 



Pfc and pk+i are density matrices of Gibbs configura- 
tions and thus of full rank. Consequently, the quantity 
S{ P k+i\\pk) + S(p k \\p k+1 ) is finite [3S]. From Eq. §B3) 
and (B6) it then follows that A(T k . n ) converges to 
AS(p k ,p k+ i) for n -t oo, i.e. 



AS(p k) p k+1 ) = lim A(7fc. n ). 



(B7) 



In other words by augmenting the intermediate points of 
which connects c k {p k ) and Cfc+i(/3fc+i) we can saturate 
the associated Clausius inequality for the fc-th step of 
the trajectory T. By repeating the same procedure for 
each of the steps of T a new trajectory emerges as the 
union of the individual sequences, 



N-l 



(B8) 



k=X 



where n is the multidimensional variable 
(ni, ri2, • • • , njv_i). Lemma [I] follows from the ad- 
ditivity of A ( 42 ) and taking the limit of each n k — » oo 
p7b. ■ 
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